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Leisingera methylohalidivorans Schaefer et al. 2002 emend. Vandecandelaere ef al. 2008 is the 
type species of the genus Leisingera. The genus belongs to the Roseobacter clade 
(Rhodobacteraceae, Alphaproteobacteria), a widely distributed lineage in marine environments. 
Leisingera and particularly L. methylohalidivorans strain MB2 T is of special interest due to its 
methylotrophy. Here we describe the complete genome sequence and annotation of this bacte- 
rium together with previously unreported aspects of its phenotype. The 4,650,996 bp long ge- 
nome with its 4,515 protein-coding and 81 RNA genes consists of three replicons, a single 
chromosome and two extrachromosomal elements with sizes of 221 kb and 285 kb. 



Introduction 



Strain MB2 T (= DSM 14336 T = ATCC BAA-92 T ) is 
the type strain of the species L. 
methylohalidivorans. L. methylohalidivorans MB2 T 
was isolated from a tide pool off the coast of Cali- 
fornia and first described in 2002 by Schaefer et 
al. [1]. The species was emended by Martens et al. 
[2] and by Vandecandelaere etal. [3]. 



versatility [5-8]. The genus name is derived in 
honor of Thomas Leisinger, on the occasion of his 
retirement and for his contributions to the under- 
standing of the biochemistry of bacterial methyl 
halide metabolism. Leisingera comprises organ- 
isms associated with their ability to grow by oxi- 
dation of methyl groups of methionine and, at 
least for L. methylohalidivorans, by oxidation of 
methyl halides as a sole energy and carbon source 
[1]. Methyl halide-degrading bacteria potentially 
play an important role in mitigating ozone deple- 
tion resulting from methyl chloride and methyl 
bromide emissions [9]. 



L. methylohalidivorans [1] is the type species of 
the genus Leisingera, which currently contains 
two more validly named species, L. aquimarina [3] 
and L. nanhaiensis [4]. The genus belongs to the 
Roseobacter clade, a widely distributed lineage in 
marine habitats with considerable metabolic 




The Genomic Standards Consortium 



Buddruhs et al. 



Here we present a summary classification and a 
set of features for L. methylohalidivorans MB2 T , 
including novel aspects of its phenotype, together 
with the description of the complete genomic se- 
quencing and annotation. 

Classification and features 
16S rRNA analysis 

Figure 1 shows the phylogenetic neighborhood of 
L. methylohalidivorans DSM 14336 T in a 16S rRNA 
based tree. The sequences of the five 16S rRNA 
gene copies in the genome differ from each other 
by up to two nucleotides, and differ by up to four 
nucleotides from the previously published 16S 
rRNA sequence (AY005463] [1]. 

A representative genomic 16S rRNA sequence of L. 
methylohalidivorans DSM 14336 T was compared 
with the Greengenes database for determining the 
weighted relative frequencies of taxa and (trun- 
cated] keywords as previously described [10]. The 
most frequently occurring genera were Ruegeria 
(32.5%], Phaeobacter (28.2%], Roseobacter 
(14.2%), Silicibacter (12.9%) and Nautella (3.5%) 
(143 hits in total). Regarding the three hits to se- 
quences from the species, the average identity 
within HSPs was 99.9%, whereas the average cov- 
erage by HSPs was 99.1%. Regarding the single hit 
to sequences from other species of the genus, the 
average identity within HSPs was 99.4%, whereas 
the average coverage by HSPs was 99.8%. Among 
all other species, the one yielding the highest 
score was 'Leisingera aquamarina' (AM900415; a 
misnomer for L. aquimarind) [3], which corre- 
sponded to an identity of 99.4% and an HSP cov- 
erage of 99.8%. (Note that the Greengenes data- 
base uses the INSDC (= EMBL/NCBI/DDBJ) anno- 
tation, which is not an authoritative source for 
nomenclature or classification.) The highest- 
scoring environmental sequence was AY007684 
('marine isolate JP88.1'), which showed an identi- 
ty of 98.1% and an HSP coverage of 100.1%. The 
most frequently occurring keywords within the 
labels of all environmental samples which yielded 
hits were 'microbi' (4.1%), 'marin' (2.8%), 
'structur' (2.3%), 'biofilm' (2.1%) and 'swro' 
(2.1%) (100 hits in total). Environmental samples 
which yielded hits of a higher score than the high- 
est scoring species were not found. This indicates 
that the species is rarely detected in the environ- 
ment. 



Morphology and physiology 

The characteristics of strain MB2 T are summa- 
rized in Table 1. Cells of L. methylohalidivorans 
MB2 T are Gram-negative and motile, obligatory 
aerobic and rod-shaped or rather pleomorphic, 
depending on the cultivation medium (Table 1) 
[1]. Colonies are non-pigmented, smooth, with an 
entire edge when grown on solid media regard- 
less of the carbon source [1]. The strain forms 
single or paired rods (1.1-1.4 x 0.4-0.5 \im) 
when grown with methyl halides, methionine or 
DMS on mineral medium. When cultured with 
yeast extract or glycine betaine, the rods become 
enlarged and elongated (2.4-8.2 x 0.7-0.8 |im). 
Yeast-grown cell lines returned to mineral salts 
medium with MeBr as the substrate reestablish 
their original form [1]. Cells grown on marine 
broth showed the standard ovoid rod morpholo- 
gy (Figure 2). 

Growth also occurs on casamino acids and weakly 
on TSA. No growth was observed on NA, R2A, PYG, 
carbon sources, amino acids (other than methio- 
nine) and small organic acids. Cells are catalase- 
and oxidase-positive. The strain does not hydro- 
lyze starch, aesculin or gelatin, and tested positive 
for leucine arylamidase activity; weak valine 
arylamidase and naphtol-AS-BI-phosphohydrolase 
activities. No activity is detected for alkaline 
phosphatase esterase (C4), esterase lipase (C8), 
lipase (C14), cystine arylamidase, trypsin, a- 
chymotrypsin, acid phosphatase, a-galactosidase, 
(B-galactosidase, (B-glucuronidase, a-glucosidase, 
iV-acetyl-(3-glucosaminidase, a-mannosidase, a- 
fucosidase, arginine dihydrolase or urease. It is 
unable to use nitrate as an electron acceptor. Vit- 
amins are not necessary for growth. Strain MB2 T 
does not degrade tyrosine, casein or DNA. No 
indole production or fermentation of glucose were 
detected [1,3]. As a marine bacterium isolated 
from seawater, growth occurred over a salinity 
range of 10-60 g/L NaCl , with an optimum at the 
salinity of seawater. The optimum Mg 2+ concentra- 
tion for strain MB2 T was 40-80 mM, which over- 
laps with the 54 mM concentration found in sea- 
water [1]. 

Strain MB2 T is susceptible to penicillin G (50 [ig), 
cefoxitin (30 [ig), erythromycin (15 ^ig), strepto- 
mycin (25 [ig) and tetracycline (30 \ig). It is mod- 
erately susceptible to gentamicin (10 \ig] but re- 
sistant to vancomycin (30 \ig], trimethoprim 
(1.25 [ig) and clindamycin (2 \ig] [1,3]. 
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- Phaeobacter daeponensis (DQ981486) 



73/69 



Phaeobacter caeruleus (AM943630) * 

Leisingera methylohalidivorans (IMG251 1206947) * 

66/85 

Leisingera aquimarina (AM90041 5) * 



Phaeobacter inhibens (AY177712) * 

Phaeobacter arcticus (DQ514304) 



Phaeobacter gallaeciensis (Y13244) * 

Leisingera nanhaiensis (FJ232451 ) * 

Seohaeicola saemankumensis (EU221 274) 

I Oceanicola nitratireducens (EU581832) 



99/99 



81/- 



78/75 



85/79 



Oceanicola nanhaiensis (DQ414420) 
Oceanicola pacificus (DQ659449) 



-Oceanicola batsensis (AY424898) 



Oceanicola granulosus (A Y424896) 



Oceanicola marinus (DQ822569) 



Figure 1. Phylogenetic tree highlighting the position of L. methylohalidivorans relative to the type strains of the other species 
within the genus Leisingera and the neighboring genera Phaeobacter and Oceanicola. The tree was inferred from 1,385 
aligned characters of the 16S rRNA gene sequence under the maximum likelihood (ML) criterion as previously described 
[10]. Oceanicola were included in the dataset for use as outgroup taxa. The branches are scaled in terms of the expected 
number of substitutions per site. Numbers adjacent to the branches are support values from 1,000 ML bootstrap replicates 
(left) and from 1,000 maximum-parsimony bootstrap replicates (right) if larger than 60% [10]. Lineages with type strain ge- 
nome sequencing projects registered in GOLD [1 1] are labeled with one asterisk, those also listed as standard 'Complete and 
Published' with two asterisks [12-14]. The genomes of three more Leisingera and Phaeobacter species are published in the 
current issue of Standards in Genomic Science [1 5-1 7]. 
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The utilization of carbon compounds by L. 
methylohalidivorans DSM 14336 T was also deter- 
mined for this study using Generation-Ill 
microplates in an OmniLog phenotyping device 
(BIOLOG Inc., Hayward, CA, USA). The microplates 
were inoculated at 28°C with a cell suspension at a 
cell density of 95-96% turbidity and dye IF-A. Fur- 
ther additives were vitamin, micronutrient and sea 
salt solutions. The exported measurement data were 
further analyzed with the opm package for R [28,29], 
using its functionality for statistically estimating pa- 
rameters from the respiration curves such as the 
maximum height, and automatically translating the- 
se values into negative, ambiguous, and positive re- 
actions. The strain was studied in two independent 
biological replicates, and reactions with a different 
behavior between the two repetitions were regard- 
ed as ambiguous and are not listed below. 

The strain gave positive reactions for 1% NaCl, 4% 
NaCl, D-glucose, D-mannitol, D-glucose-6- 
phosphate, D-aspartic acid, L-alanine, L-arginine, L- 
glutamic acid, L-histidine, L-pyroglutamic acid, L- 
serine, D-galacturonic acid, glucuronamide, quinic 
acid, D-saccharic acid, D-lactic acid methyl ester, a- 
keto-glutaric acid, L-malic acid, nalidixic acid, 
acetoacetic acid, propionic acid and acetic acid. 

The strain was negative for sucrose, pH 6, pH 5, D- 
melibiose, D-salicin, N-acetyl-D-glucosamine, N- 
acetyl-D-galactosamine, N-acetyl-neuraminic acid, 
8% NaCl, D-galactose, 3-0-methyl-D-glucose, D- 
fucose, L-fucose, inosine, 1% sodium lactate, fusidic 
acid, D-serine, D-sorbitol, D-arabitol, D-fructose-6- 
phosphate, D-serine, troleandomycin, rifamycin SV, 
minocycline, lincomycin, guanidine hydrochloride, 
niaproof 4, pectin, L-galactonic acid-y-lactone, mucic 
acid, vancomycin, tetrazolium violet, tetrazolium 
blue, p-hydroxy-phenylacetic acid, methyl pyruvate, 
citric acid, bromo-succinic acid, potassium tellurite, 
a-hydroxy-butyric acid, (B-hydroxy-butyric acid, a- 
keto-butyric acid, sodium formate, aztreonam, bu- 
tyric acid and sodium bromate. 

Regarding the common subset of growth experi- 
ments and OmniLog experiments, the results were 
identical with few exceptions. Expectedly [30], on 
some substrates respiration was detected by phe- 
notype microarray analysis even though these 
substrates did not sustain growth. 

Chemotaxonomy 

The major respiratory lipoquinone present is Q10 
[1]. The polar lipids comprise phosphatidylglycerol, 



phosphatidylethanolamine, an unidentified phos- 
pholipid, two unidentified lipids and an aminolipid. 
Phosphatidylcholine is not present. The fatty acids 
comprise Cio : o 3-oh, Cu-.i, C12.0 3-oh, Ci6:o, Ci6 : o 2-oh, 

Cl8:la)9c, Cl8:lco7c, Cl8:0 and 11-methyl C 18:10)70 The Cl0:0 

3-oh and Ci6:o 3-oh fatty acids are ester-linked, while 
the Ci2:0 3-oh fatty acid is amide-linked [2]. 

Genome sequencing and annotation 

Genome project history 

This organism was selected for sequencing on the 
basis of the DOE Joint Genome Institute Community 
Sequencing Program (CSP) 2010, CSP 441 "Whole 
genome type strain sequences of the genera 
Phaeobacter and Leisingera - a monophyletic group 
of physiological highly diverse organisms". The ge- 
nome project is deposited in the Genomes On Line 
Database [11] and the complete genome sequence 
is deposited in GenBank and the Integrated Micro- 
bial Genomes database (IMG] [31]. Sequencing, fin- 
ishing and annotation were performed by the DOE 
Joint Genome Institute QGI] using state of the art 
sequencing technology [32]. A summary of the pro- 
ject information is shown in Table 2. 

Growth conditions and DNA isolation 

A culture of DSM 14336 T was grown aerobically in 
DSMZ medium 514 [33] at 20°C. Genomic DNA was 
isolated using a Jetflex Genomic DNA Purification 
Kit (GENOMED 600100] following the standard 
protocol provided by the manufacturer but modi- 
fied by an incubation time of 40 min, the incubation 
on ice over night on a shaker, the use of additional 
10 [il proteinase K, and the addition of 100 [il pro- 
tein precipitation buffer. DNA is available from 
DSMZ through the DNA Bank Network [34]. 

Genome sequencing and assembly 

The draft genome sequence was generated using 
Illumina sequencing technology. For this genome, 
we constructed and sequenced an Illumina short- 
insert paired-end library with an average insert size 
of 270 bp, which generated 10,989,662 reads. In ad- 
dition, an Illumina long-insert paired-end library 
with an average insert size of 9,000 bp was con- 
structed, generating 1,005,012 reads for a total of 
1,798 Mb of Illumina data (Feng Chen, unpublished]. 
All general aspects of library construction and se- 
quencing performed can be found at the JGI web site 
[35]. The initial draft assembly contained 16 contigs 
in 6 scaffoldfs]. The initial draft data was assembled 
with Allpaths [36] and the consensus was computa- 
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tionally shredded into 10 kbp overlapping fake 
reads (shreds]. The Illumina draft data was also as- 
sembled with Velvet [37], and the consensus se- 
quences were computationally shredded into 1.5 
kbp overlapping fake reads (shreds]. The Illumina 
draft data was assembled again with Velvet using 
the shreds from the first Velvet assembly to guide 
the next assembly. The consensus from the second 
Velvet assembly was shredded into 1.5 kbp overlap- 
ping fake reads. The fake reads from the Allpaths 
assembly and both Velvet assemblies and a subset of 
the Illumina CLIP paired-end reads were assembled 



using parallel phrap (High Performance Software, 
LLC] [38]. Possible mis-assemblies were corrected 
with manual editing in Consed [38]. Gap closure was 
accomplished using repeat resolution software (Wei 
Gu, unpublished], and sequencing of bridging PCR 
fragments with Sanger technologies. A total of 15 
additional sequencing reactions were completed to 
close gaps and to raise the quality of the final se- 
quence. The total size of the genome is 4,630,996 bp 
and the final assembly is based on 1,798 Mb of 
Illumina draft data, which provides an average 382.5 
x coverage of the genome. 



Table 1. Classification and general features of L. methylohalidivorans MB2 T according to the MIGS recommenda- 
tions [18] published by the Genome Standards Consortium [19]. 



MIGS ID Property 



Term 



Evidence code 







Domain Bacteria 


TAS [20] 






Phylum Proteobacteria 


TAS [21] 






Class Alphaproteobacteria 


TAS [22,23] 




Current classification 


Order Rhodobacterales 
Family Rhodobacteraceae 
Genus Leisingera 

Species Leisingera methylohalidivorans 


TAS [23,24] 
TAS [23,25] 
TAS [1-3] 
TAS [1,3] 


MIGS-7 


Subspecific genetic lineage (strain) 


MB2 T 


TAS [1] 


MIGS-12 


Reference for biomaterial 


Schaefer et al. 2002 


TAS [1] 




Gram stain 


negative 


TAS [1] 




Cell shape 


ovoid rods/ pleomorphism 


TAS [1] 




Motility 


motile 


TAS [1] 




Sporulation 


non-sporulating 


TAS [1] 




Temperature range 


mesophile 


TAS [1] 




Optimum temperature 


27°C 


TAS [1] 




Salinity 


halophile 


TAS [1] 


MIGS-22 


Relationship to oxygen 


obligatory aerobic 


TAS [1] 




Carbon source 


complex substrates, methyl halides, DMS, 
methionine, glycine betaine 


TAS [1] 


MIGS-6 


Habitat 


aquatic, sea water 


TAS [1] 


MIGS-6.2 


P H 


7.7 


TAS [1] 


MIGS-15 


Biotic relationship 


free-living 


TAS [1] 


MIGS-14 


Known pathogenicity 


none 


TAS [1] 




Biosafety level 


1 


TAS [26] 


MIGS-23.1 


Isolation 


tide pool 


TAS [1] 


MIGS-4 


Geographic location 


USA, Washington DC 


TAS [1] 


MIGS-5 


Time of sample collection 


2002 or before 


TAS [1] 


MIGS-4. 1 


Latitude 


38.90 


TAS [1] 


MIGS-4.2 


Longitude 


-77.03 


TAS [1] 



Evidence codes - IDA: Inferred from Direct Assay; TAS: Traceable Author Statement (i.e., a direct report exists in the lit- 
erature); NAS: Non-traceable Author Statement (i.e., not directly observed for the living, isolated sample, but based on a 
generally accepted property for the species, or anecdotal evidence). Evidence codes are from the Gene Ontology project 
[27]. 
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Table 2. Genome sequencing project information 



MIGS ID 


Property 


Term 


M1GS-31 


Finishing quality 


Finished 


MIGS-28 


Libraries used 


Two lllumina paired-end libraries (270 bp and 9 kb insert size) 


MIGS-29 


Sequencing platforms 


lllumina GAii 


MIGS-31.2 


Sequencing coverage 


382.5 x lllumina 


MIGS-30 


Assemblers 


Allpaths, Velvet 1 .1 .05, phrap version SPS - 4.24 


Mir^Q 1.1 

/Vllvj j- 3Z 


Gene calling method 


rrOulgdl I .4, vjGnerKI/Vlr 




INSDC ID 


INSDC ID CP006773 (cMeth_4145), CP006774 (pMeth_B221 ), 
CP006775 (pMeth_A285) 




GenBank Date of Release 


September 30, 201 3 




GOLD ID 


Gi10858 






PP INI A7A^71 




Database: IMG-GEBA 


2512564009 


MIGS-13 


Source material identifier 


DSM 14336 




Project relevance 


Tree of Life, carbon cycle, sulfur cycle, environmental 



Genome annotation 

Genes were identified using Prodigal [39] as part 
of the DOE-JGI genome annotation pipeline [40], 
followed by a round of manual curation using the 
JGI GenePRIMP pipeline [41]. The predicted CDSs 
were translated and used to search the National 
Center for Biotechnology Information (NCBI] 
nonredundant database, UniProt, TIGR-Fam, Pfam, 
PRIAM, KEGG, COG, and InterPro databases. Addi- 
tional gene prediction analysis and functional an- 
notation was performed within the Integrated Mi- 
crobial Genomes - Expert Review (IMG-ER] plat- 
form [31]. 

Genome properties 

The L. methylohalidivorans DSM 14336 T genome 
statistics are provided in Table 3 and Figure 3. The 
genome consists of three circular replicons with a 
total length of 4,650,996 bp and a G+C content of 
62.3%. The replicons correspond to a single 
chromosome (4,144,900 bp in length] and two 
extrachromosomal elements 220,701 bp and 
285,395 bp in length. Of the 4,596 genes predict- 
ed, 4,515 were protein-coding genes, and 81 
RNAs. In addition, 293 pseudogenes were also 



identified. The majority of the protein-coding 
genes (77.4%] were assigned a putative function 
while the remaining ones were annotated as hypo- 
thetical proteins. The distribution of genes into 
COGs functional categories is presented in Table 4. 

Insights into the genome 

Genome sequencing of L. methylohalidivorans DSM 
14336 T reveals the presence of two plasmids with 
sizes of 221 kb and 285 kb (Table 5]. The circular 
conformation of the chromosome and the two 
extrachromosomal elements have been experi- 
mentally validated. The two plasmids contain 
characteristic replication modules of the DnaA- 
like and RepABC-type comprising a replicase as 
well as the parAB partitioning operon [42]. The 
respective replicases that mediate the initiation of 
replication are designated according to the estab- 
lished plasmid classification scheme [43]. The dif- 
ferent numbering of the replicase RepC-8 from the 
RepABC-type plasmids corresponds to specific 
plasmid compatibility groups that are required for 
a stable coexistence of the replicons within the 
same cell [44]. 
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Figure 3. Graphical map of the chromosome (cMeth_4145) and the two extrachromosomal elements 
(pMeth_B221 and 75.2). From outside to the center: Genes on forward strand (color by COG catego- 
ries), Genes on reverse strand (color by COG categories), RNA genes (tRNAs green, rRNAs red, other 
RNAs black), GC content, GC skew. 
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Table 3. Genome statistics 



Attribute 


Number 


% of Total 


Genome size (bp) 


4,650,996 


100.00 


DNA coding region (bp) 


3,929,972 


84.50 


DNA G+C content (bp) 


2,898,874 


62.33 


Number of replicons 


3 




Extrachromosomal elements 


2 




Total genes 


4,596 


100.00 


RNA genes 


81 


1.76 


rRNA operons 


5 




tRNA genes 


62 


1.35 


Protein-coding genes 


4,515 


98.24 


Pseudo genes 


293 


6.38 


Genes with function prediction 


3,558 


77.42 


Genes in paralog clusters 


1,675 


36.44 


Genes assigned to COGs 


3,493 


76.00 


Genes assigned Pfam domains 


3,567 


77.61 


Genes with signal peptides 


1,470 


31.98 


Genes with transmembrane helices 


867 


18.86 


CRISPR repeats 


0 





Table 4. Number of genes associated with the general COG functional categories 



Code 


Value 


%age 


Description 


J 


179 


4.67 


Translation, ribosomal structure and biogenesis 


A 


0 


0 


RNA processing and modification 


K 


322 


8.39 


Transcription 


L 


232 


6.05 


Replication, recombination and repair 


B 


3 


0.08 


Chromatin structure and dynamics 


D 


34 


0.89 


Cell cycle control, cell division, chromosome partitioning 


Y 


0 


0 


Nuclear structure 


V 


40 


1.04 


Defense mechanisms 


T 


161 


4.2 


Signal transduction mechanisms 


M 


181 


4.72 


Cell wall/membrane/envelope biogenesis 


N 


51 


1.33 


Cell motility 


Z 


1 


0.03 


Cytoskeleton 


W 


0 


0 


Extracellular structures 


U 


62 


1.62 


Intracellular trafficking, secretion, and vesicular transport 


o 


138 


3.6 


Posttranslational modification, protein turnover, chaperones 


c 


237 


6.18 


Energy production and conversion 


G 


154 


4.01 


Carbohydrate transport and metabolism 


E 


459 


11.96 


Amino acid transport and metabolism 


F 


102 


2.66 


Nucleotide transport and metabolism 


H 


184 


4.8 


Coenzyme transport and metabolism 


I 


150 


3.91 


Lipid transport and metabolism 


P 


182 


4.74 


Inorganic ion transport and metabolism 


Q 


123 


3.21 


Secondary metabolites biosynthesis, transport and catabolism 


R 


465 


12.12 


General function prediction only 


S 


377 


9.83 


Function unknown 




1110 


24.09 


Not in COGs 
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Table 5. General genomic features of the chromosome and extrachromosomal replicons from L. 
methylohalidivorans strain DSM 14336 T . 



Replicon 



Scaffold Replicase Length (bp) GC (%) Topology 



No. Genes 



cMeth_4145 


1 


DnaA 


4,144,900 


62 


circular 


4,135 


pMeth_A285 


2 


DnaA-like I 


285,395 


62 


circular 


269 


pMeth_B221 


3 


RepC-8 


220,701 


62 


circular 


204 



+ deduced from automatic annotation. 



Table 6. Integrated Microbial Genome (IMG) locus tags of L. methylohalidivorans DSM 14336 T genes + 
Replicon Replication initiation Plasmid stability Type IV secretion 



cMeth_4145 


Replicase 

DnaA 


Locus Tag 

Meth_0476 


Toxin 


Antitoxin 


VirB4 VirD4 


pMeth_A285 


DnaA-like I 


Meth_0245 


Meth_0528 


Meth_0529 




pMeth_B221 


RepC-8 


Meth_0107 






Meth_0035 1 



+ Genes for the initiation of replication, toxin/antitoxin modules and type IV secretion systems (T4SS) that are 
required for conjugation. 

1 Presence of adjacent DNA relaxase VirD2. 



The locus tags of all replicases, plasmid stability 
modules and the large virB4 and virD4 genes of 
the type IV secretion systems are presented in Ta- 
ble 6. The larger plasmid, pMeth_A285, harbors a 
postsegregational killing system (PSK] consisting 
of a typical operon with two small genes encoding 
a stable toxin and an unstable antitoxin [45]. The 
smaller plasmid pMeth_B221 contains the virD2 
and virD4 genes of the type IV secretion system, 
but it is probably non-conjugative, since the virB 
operon for the formation of a transmembrane 
channel is missing [46,47]. 

The 285 kb DnaA-like I replicon pMeth_A285 con- 
tains a large type VI secretion system (T6SS] with 
a size of about 30 kb. The role of this export sys- 
tem was first described in the context of bacterial 
pathogenesis, but recent findings indicate a more 
general physiological role in defense against eu- 
karyotic cells and other bacteria in the environ- 
ment [48]. Homologous T6S systems are present 
on the DnaA-like I plasmids of L. aquimarina DSM 
24565 T (pAqui_F126] and Phaeobacter caeruleus 
DSM 24564? (pCaer_C109] as well as the RepC-8 
type plasmid of Phaeobacter daeponensis DSM 
23529 T (pDaep_A276]. This extrachromosomal 



replicon also harbors a TonB-dependent 
siderophore receptor (Meth_0471] and genes of a 
putative ABC-type Fe 3+ siderophore transport sys- 
tem (Meth_0472 to Meth_0467]. 

The 221 kb RepC-8 type replicon pMeth_B221 
contains five ABC-transporters. One of them, 
which probably transports nitrate/sulfonate or 
bicarbonate (Meth_0002, Meth_0001, Meth_0204, 
Meth_0203], is located adjacent to the large and 
small subunit genes of the nitrate reductase (EC 
1.7.1.4; Meth_0202, Meth_0201] and an anaerobic 
dehydrogenase (EC 1.7.99.4; Meth_0200] hence 
indicating a functional role of the plasmid in an- 
aerobic metabolism. 

To quantify the differences in COG functional cate- 
gories between the three replicons and to deter- 
mine the over-represented categories, we used 
approaches based on entropy and the broken- 
stick distribution, respectively. We applied these 
methods to all genes that were assigned to a COGs 
category from either genome [49]. Figure 4 shows 
the bar plot of the COG categories of the replicons 
[18]. The analysis revealed one over-represented 
COG category for the small extrachromosomal el- 
ement (pMeth_B221], i.e. "amino acid metabolism" 
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(category E]. For instance, this replicon encodes 
nine spermidine/putrescine transporter sequenc- 
es (Meth_0060, _0061, _0062, _0063, _0133, _0134, 
_0135, _0136, _0169] suggesting that these com- 
pounds are an important source for L. 
methylohalidivorans DSM 14336 T . Spermidine and 
putrescine are produced in marine phytoplankton 
and zooplankton to regulate cell proliferation and 
bloom formation [50]. 

The COG category P ("inorganic ion transport and 
metabolism"] (Figure 4) is highly represented in 



the larger extrachromosomal element 
(Meth_0238, _0261, _0263, _0264, _0265, _0266, 
_0303, _0305, _0355, _0360, _0378, _0413, _0414, 
_0415, _0463, _0468, _0469, _0470, _0471). This 
replicon encodes a broad spectrum of inorganic 
transport and regulation systems for sulfate, 
phosphate, 2-aminoethylphosphate, manga- 
nese(II], zinc(II], ferric, ferrous, ferric-citrate, 
formate, nitrite, calcium(II], sodium, molybdenum 
and copper. 
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Figure 4. Bar plot of the relative amounts of the COG categories of the chromosome (cMeth_4145 = 
Leime_Contig76.3, left) and both extrachromosomal elements (pMeth_A285 = Leime_Contig 75.2, center, and 
pMeth_ B221 = Leime_Contig74.1, right) The COG functional categories are described in Table 4. 



In accordance with the known ability of L. 
methylohalidivorans DSM 14336 T to grow by oxi- 
dation of methyl halides [1], the genome analysis 
revealed the genes for the proposed pathway of 
methyl chloride metabolism as described by 
McDonald et al. 2002 [9]. Using the JGI-IMG 
BLASTp tool [51,52], the gene for first 
methyltransferase I (cmuA) indeed yielded a hit to 
the gene cmuA ("predicted cobalamin binding pro- 
tein", Meth_2531] in the genome of L. 



methylohalidivorans DSM 14336 T , with a sequence 
similarity of 31%. Searching for the second en- 
zyme methyltransferase II [cmuE) yielded a hit to 
the enzyme adjacent to the predicted cobalamin 
binding protein ("methionine synthase I 
(cobalamin-dependent), methyltransferase do- 
main", Meth_2528]. For the next enzymes in the 
methyl-chloride metabolism, we compared the 
genes metF,folD, purU and FDH and found the fol- 
lowing results: 39% similarity to a 5,10- 
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methylenetetrahydrofolate reductase 
(Meth_1763), for metF; 56% to a 5,10-methylene- 
tetrahydrofolate dehydrogenase/Methenyl- 
tetrahydrofolate cyclohydrolase (Meth_4077, 
Meth_3180] for folD; 36% to a 
phosphoribosylglycinamide formyltransferase, 
formyltetrahydrofolate-dependent (Meth_2536] 
for purU; and 79% to a formate dehydrogenase 
(Meth_4011] for FDH. 

An estimate of the DNA-DNA hybridization (DDH] 
similarity between L. methylohalidivorans DSM 
14336 T and the draft genomes of L. aquimarina 
DSM 24565 T , L. nanhaiensis DSM 24252 T , P. 
arcticus DSM 23566 T , P. caeruleus DSM 24564 T , P. 
daeponensis DSM 23529 T , P. gallaeciensis CIP 
105210 T and P. inhibens DSM 16374 T was gener- 
ated with the GGDC Genome-to-Genome Distance 
Calculator version 2.0 [53-55]. This system calcu- 
lates the distances by comparing the genomes to 
obtain HSPs (high-scoring segment pairs] and 



interfering distances from three formulae (1, HSP 
length / total length; 2, identities / HSP length; 3, 
identities / total length] [53]. Table 7 shows the 
results of the pairwise comparisons between L. 
methylohalidivorans DSM 14336 T and the other 
seven genomes. As the results of the 16S rRNA 
analysis (Figure 1] revealed, the two Leisingera 
species L. methylohalidivorans and L. aquimarina 
show a close relationship, whereas L. nanhaiensis 
does not cluster together with the other two 
Leisingera species. The DDH similarities calculated 
in silico yielded similar results, indicating that the 
classification of L. nanhaiensis might need to be 
reconsidered. Furthermore, the DDH similarities 
of L. methylohalidivorans to Phaeobacter species 
are not significantly smaller, especially in the case 
of P. caeruleus and P. daeponensis, than to L. 
aquimarina and as already described to L. 
nanhaiensis. 



Table 7. DDH similarities between L methylohalidivorans DSM 14336 T and the other Leisingera and Phaeobacter spe- 



cies (with genome-sequenced type strains) calculated in 


silico with the GGDC server version 2.0 [55]. 


Reference species 


HSP length / 
total length [%] 


identities / HSP length [%] 


identities / total length [%] 


L. aquimarina (AXBE00000000) 


52.40 ± 3.47 


32.40 ± 2.46 


47.00 ± 3.03 


L nanhaiensis (AXBG00000000) 


14.50 ± 3.11 


19.20 ± 2.29 


14.60 ± 2.64 


P. arcticus (AXBF00000000) 


17.20 ± 3.28 


20.40 ± 2.32 


17.00 ± 2.77 


P. caeruleus (AXBI00000000) 


45.80 ± 3.41 


27.00 ± 2.42 


39.90 ± 3.01 


P. daeponensis (AXBD00000000) 


48.70 ± 3.43 


26.90 ± 2.42 


41.90 ± 3.01 


P. gallaeciensis (AOQA01 000000) 


17.90 ± 3.31 


21.00 ± 2.33 


17.60 ± 2.80 


P. inhibens (AXBB00000000) 


18.50 ± 3.34 


21.10 ± 2.33 


18.10 ± 2.82 



The standard deviations indicate the inherent uncertainty in estimating DDH values from intergenomic distances 
based on models derived from empirical test data sets (which are always limited in size); see [55] for details. The dis- 
tance formulae are explained in [54]. The numbers in parentheses are GenBank accession numbers identifying the 
underlying genome sequences. 



Acknowledgements 

The authors would like to gratefully acknowledge the 
assistance of Iljana Schroeder for growing L. 
methylohalidivorans cultures and Evelyne-Marie 
Brambilla for DNA extraction and quality control (both 
at the DSMZ]. The work conducted by the U.S. Depart- 
ment of Energy Joint Genome Institute was supported 
by the Office of Science of the U.S. Department of Ener- 
gy under contract No. DE-AC02-05CH11231; AL was 
supported by Russian Ministry of Science Mega-grant 



no.ll.G34.31.0068; SJ O'Brien Principal Investigator. 
The work conducted by the members of the 
Roseobacter consortium was supported by the German 
Research Foundation (DFG] Transregio-SFB 51 with 
PhD stipends for NB and AF. We also thank the Europe- 
an Commission which supported phenotyping via the 
Microme project 222886 within the Framework 7 pro- 
gram. 



138 



Standards in Genomic Sciences 



Buddruhs et al. 



References 

1. SchaeferJK, Goodwin KD, Mcdonald IR, Murrell 
JC, Oremland RS. Leisingera methylohalidivorans 
gen. nov., sp. nov., a marine methylotroph that 
grows on methyl bromide. Int J Syst Evol Microbiol 
2002; 52:851-859. PubMed 
http://dx.doi.Org/1 0.1 099/ijs.0.01 960-0 

2. Martens T, Heidorn T, Pukall R, Simon M, Tindall 
BJ, Brinkhoff T. Reclassification of Roseobacter 
gallaeciensis Ruiz-Ponte et al. 1 998 as Phaeobacter 
gallaeciensis gen. nov., comb, nov., description of 
Phaeobacter inhibens sp. nov., reclassification of 
Ruegeria algicola (Lafay ef al. 1 995) Uchino et al. 

1 999 as Marinovum algicola gen. nov., comb, 
nov., and emended descriptions of the genera 
Roseobacter, Ruegeria and Leisingera. Int J Syst 
Evol Microbiol 2006; 56:1 293-1 304. PubMed 
http://dx.doi.Org/1 0.1 099/ijs.0.63724-0 

3. Vandecandelaere I, Segaert E, Mollica A, Faimali 
M, Vandamme P. Leisingera aquimarina sp. nov., 
isolated from a marine electroactive biofilm, and 
emended descriptions of Leisingera 
methylohalidivorans Schaefer et al. 2002, 
Phaeobacter daeponensis Yoon et al. 2007 and 
Phaeobacter inhibens Martens et al. 2006. Int J Syst 
Evol Microbiol 2008; 58:2788-2793. PubMed 
http://dx.doi.Org/1 0.1 099/ijs.0.65844-0 

4. Sun F, Wang B, Liu X, Lai Q, Du Y, Li G, Luo J, 
Shao Z. Leisingera nanhaiensis sp. nov., isolated 
from marine sediment. Int J Syst Evol Microbiol 
2010; 60:275-280. PubMed 
http://dx.doi.Org/1 0.1 099/ijs.0.01 0439-0 

5. Berger M, Brock NL, Liesegang H, Dogs M, Preuth 
I, Simon M, DickschatJS, Brinkhoff T. Genetic 
analysis of the upper phenylacetate catabolic 
pathway in the production of tropodithietic acid by 
Phaeobacter gallaeciensis. Appl Environ Microbiol 
2012; 78:3539-3551 . PubMed 
http://dx.doi.Org/10.1 128/AEM.07657-1 1 

6. Brinkhoff T, Giebel HA, Simon M. Diversity, ecol- 
ogy, and genomics of the Roseobacter clade: a 
short overview. Arch Microbiol 2008; 189:531 - 
539. PubMed http://dx.doi.org/10.1007/s00203- 
008-0353-v 

7. Slightom RN, Buchan A. MINIREVIEW Surface 
colonization by marine Roseobacters: Integrating 
genotype and phenotype. Appl Environ Microbiol 
2009; 75:6027-6037. PubMed 
http://dx.doi.Org/1 0.1 1 28/AEM.01 508-09 

8. Wagner- Dob ler I, Biebl H. Environmental biology 
of the marine Roseobacter lineage. Annu Rev 
Microbiol 2006; 60:255-280. PubMed 



http://standardsingenomics.org 



http://dx.doi.Org/10.1 146/annurev. micro. 60. 08080 
5.142115 

9. McDonald IR, Warner KL, Mcanulla C, Woodall 
CA, Oremland RS, Murrell JC. A review of bacteri- 
al methyl halide degradation: biochemistry, genet- 
ics and molecular ecology. Environ Microbiol 
2002; 4:193-203. PubMed 

http ://dx .doi .org/1 0 . 1 046/j . 1 462 - 
2920.2002.00290.x 

1 0. Goker M, Cleland D, Saunders E, Lapidus A, Nolan 
M, Lucas S, Hammon N, Deshpande S, Cheng JF, 
Tapia R, et al. Complete genome sequence of 
Isosphaera pallida type strain (IS1 B T ). Stand Ge- 
nomic Sci 201 1 ; 4:63-71 . PubMed 
http://dx.doi.Org/1 0.4056/sigs.1 533840 

1 1 . Pagani I, Liolios K, Jansson J, Chen IM, Smirnova T, 
Nosrat B, Markowitz VM, Kyrpides NC. The Ge- 
nomes OnLine Database (GOLD) v. 4: status of ge- 
nomic and metagenomic projects and their associ- 
ated metadata. Nucleic Acids Res 2012; 40:D571- 
D579. PubMed 

http://dx.doi.Org/1 0.1 093/nar/gkr1 1 00 

12. Beyersmann PG, Chertkov O, Petersen J, Fiebig A, 
Chen A, Pati A, Ivanova N, Lapidus A, Goodwin 
LA, Chain P, etal. Genome sequence of 
Phaeobacter caeruleus type strain (DSM 24564 T ), a 
surface-associated member of the marine 
Roseobacter clade. Stand Genomic Sci 201 3; 
8:403-419. http://dx.doi.org/10.4056/sigs.3927626 

13. Freese H, Dalingault H, Petersen J, Pradella S, 
Fiebig A, Davenport K, Teshima H, Chen A, Pati A, 
Ivanova N, et al. Genome sequence of the plasmid 
and phage-gene rich marine Phaeobacter arcticus 
type strain (DSM 23566 T ). Stand Genomic Sci 
2013; 8:450-464. 

http://dx.doi.org/10.4056/sigs.383362 

14. Riedel T, Teshima H, Petersen J, Fiebig A, Daven- 
port K, Dalingault H, Erkkila T, Gu W, Munk C, Xu 
Y, et al. Genome sequence of the Leisingera 
aquimarina type strain (DSM 24565 T ), a member of 
the Roseobacter clade rich in extrachromosomal 
elements. Stand Genomic Sci 201 3; 8:389-402. 
http://dx.doi.org/10.4056/sigs.3858183 

15. Dogs M, Teshima H, Petersen J, Fiebig A, Chertkov 
O, Dalingault H, Chen A, Pati A, Goodwin LA, 
Chain P, et al. Genome sequence of Phaeobacter 
daeponensis strain DSM 24529 T , a facultatively 
anaerobic member of the genus Phaeobacter iso- 
lated from marine sediment. Stand Genomic Sci 
2013; 8:142-159. 

http://dx.doi.org/10.4056/sigs.4287962 



139 



Leisingera methylohalidivorans type strain (DSM 14336T) 

1 6. Dogs M, Voget S, Teshima H, Petersen J, Fiebig A, 
Davenport K, Dalingault H, Chen A, Pati A, 
Ivanova N, et al. Genome sequence of 
Phaeobacter inhibens strain T5 T , a 
secondarymetabolite producing member of the 
marine Roseobacter clade. Stand Genomic Sci 
201 3; (In press). 

17. Breider S, Teshima H, Petersen J, Fiebig A, 
Chertkov O, Dalingault H, Chen A, Pati A, Ivanova 
N, Lapidus A, et al. Genome sequence of 
Leisingera nanhaiensis strain DSM 23252 T isolated 
from marine sediment. Stand Genomic Sci 2013; 
(In press). 

18. Field D, Garrity G, Gray T, Morrison N, Selengut J, 
Sterk P, Tatusova T, Thomson N, Allen MJ, 
Angiuoli SV, etal. The minimum information 
about a genome sequence (MIGS) specification. 
Nat Biotechnol 2008; 26:541 -547. PubMed 
http://dx.doi.Org/1 0.1 038/nbt1 360 

1 9. Field D, Amaral-Zettler L, Cochrane G, Cole JR, 
Dawyndt P, Garrity GM, Gilbert J, Glockner FO, 
Hirschman L, Karsch-Mzrachi I, etal. The Ge- 
nomic Standards Consortium (GSC). PLoS Biol 
201 1 ; 9:e1 001 088. PubMed 
http://dx.doi.Org/1 0.1 371 /journal.pbio.1 001 088 

20. Woese CR, Kandler O, Wheelis ML. Towards a 
natural system of organisms: proposal for the do- 
mains Archaea, Bacteria, and Eucarya. Proc Natl 
Acad Sci USA 1990; 87:4576-4579. PubMed 
http://dx.doi.Org/1 0.1 073/pnas.87.1 2.4576 

21 . Garrity G, Bell J, Lilburn T. Phylum XIV. 
Proteobacteria phyl. nov. In: Brenner D, Krieg N, 
Staley J, Garrity G, eds. Bergey's Manual of Sys- 
tematic Bacteriology, Vol. 2 (The Proteobacteria), 
Part B (The Gammaproteobacteria). Second Edi- 
tion. New York: Springer; 2005:1 . 

22. Garrity G, Bell J, Lilburn T. Class I. 
Alphaproteobacteria class, nov. In: Garrity G, 
Brenner D, Krieg N, Staley J, eds. Bergey's Manual 
of Systematic Bacteriology, Volume 2, Part C. Se- 
cond Edition. New York: Springer; 2005:1 . 

23. Validation List No. 1 07. List of new names and 
new combinations previously effectively, but not 
validly, published. Int J Syst Evol Microbiol 2006; 
56:1-6. PubMed 

http://dx.doi.Org/1 0.1 099/ijs.0.641 88-0 

24. Garrity G, Bellm J, Lilburn T. Order III. 
Rhodobacterales ord. nov. In: Garrity G, Brenner 
D, Krieg N, Staley J, eds. Bergey's Manual of Sys- 
tematic Bacteriology, Volume 2, Part C. Second 
Edition. New York: Springer; 2005:161. 



25. Garrity GM, Bell JA, Lilburn T. Family III. 
Rhodobacteraceae fam. nov. In: Garrity GM, 
Brenner DJ, Krieg NR, Staley JT, eds. Bergey's 
Manual of Systematic Bacteriology, Volume 2, Part 
C. Second Edition. New York: Springer; 2005:161. 

26. BAuA. Classification of Bacteria and Archaea in 
risk groups. TRBA 201 0; 466:93. 

27. Ashburner M, Ball CA, Blake JA, Botstein D, Butler 
H, Cherry JM, Davis AP, Dolinski K, DwightSS, 
Eppig JT, ef al. Gene ontology: tool for the unifica- 
tion of biology. The Gene Ontology Consortium. 
Nat Genet 2000; 25:25-29. PubMed 
http://dx.doi.org/10.1038/75556 

28. Vaas LAI, Sikorski J, Hofner B, Buddruhs N, Fiebig 
A, Klenk HP, Goker M. opm: An R package for an- 
alysing OmniLog® Phenotype MicroArray Data. 
Bioinformatics 201 3; 29:1 823-1 824. PubMed 
http://dx.doi.org/10.1093/bioinformatics/btt291 

29. Vaas LAI, Sikorski J, Michael V, Goker M, Klenk 
HP. Visualization and curve-parameter estimation 
strategies for efficient exploration of phenotype mi- 
croarray kinetics. PLoS ONE 2012; 7:e34846. 
PubMed 

http://dx.doi.org/10.1371/journal.pone.0034846 

30. Vaas LAI, Marheine M, Sikorski J, Goker M, Schu- 
macher M. Impacts of pr-1 0a overexpression at the 
molecular and the phenotypic level. Int J Mol Sci 
2013; 14:15141-15166. PubMed 
http://dx.doi.Org/1 0.3390/ijms1 4071 5141 

31 . Markowitz VM, Mavromatis K, Ivanova NN, Chen 
IA, Chu K, Kyrpides NC. IMG ER: a system for mi- 
crobial genome annotation expert review and 
curation. Bioinformatics 2009; 25:2271-2278. 
PubMed 

http://dx.doi.Org/1 0.1 093/bioinformatics/btp393 

32. Mavromatis K, Land ML, Brettin TS, Quest DJ, 
Copeland A, Clum A, Goodwin L, Woyke T, 
Lapidus A, Klenk HP, et al. The fast changing land- 
scape of sequencing technologies and their impact 
on microbial genome assemblies and annotation. 
PLoS ONE 2012; 7:e48837. PubMed 
http://dx.doi.org/10.1371/journal.pone.0048837 

33. List of growth media used at the DSMZ: 
http://www.dmsz.de/catalogues/cataloque- 
microorganisms/culture-technology/list-of-media- 
for-microorganisms.html. 

34. Gemeinholzer B, Droge G, Zetzsche H, 
Haszprunar G, Klenk HP, Guntsch A, Berendsohn 
JW. The DNA Bank Network: the start from a 
German initiative. Biopres Biobanking 2011; 9:51- 
55. http://dx.doi.Org/1 0.1 089/bio.201 0.0029 



140 



Standards in Genomic Sciences 



Buddruhs et al. 



35. The DOE Joint Genome Institute. 
http ://www . j gi . doe . gov . 

36. Butler J, MacCallum I, Kleber M, Shlyakhter IA, 
Belmonte MK, Lander ES, Nusbaum CJD. 
ALLPATHS: de novo assembly of whole-genome 
shotgun microreads. Genome Res 2008; 18:810- 
820. PubMed 

http://dx.doi.org/10.1101/gr.7337908 

37. Zerbino DR, Birney E. Velvet: algorithms for de 
novo short read assembly using de Bruijn graphs. 
Genome Res 2008; 18:821-829. PubMed 
http://dx.doi.Org/1 0.1 1 01/gr.074492.1 07 

38. Phrap and Phred for Windows. MacOS, Linux, and 
Unix, http://www.phrap.com. 

39. Hyatt D, Chen G, Locascio PF, Land ML, Larimer 
FW, Hauser LJ. Prodigal: Prokaryotic gene recogni- 
tion and translation initiation site identification. 
BMC Bioinformatics 2010; 11:1 19. PubMed 
http://dx.doi.Org/1 0.1 1 86/1 471 -21 05-1 1-119 

40. Mavromatis K, Ivanova NN, Chen IM, Szeto E, 
Markowitz VM, Kyrpides NC. The DOE-JGI Stand- 
ard operating procedure for the annotations of mi- 
crobial genomes. Stand Genomic Sci 2009; 1:63- 
67. PubMed http://dx.doi.org/10.4056/sigs.632 

41 . Pati A, Ivanova NN, Mikhailova N, Ovchinnikova 
G, Hooper SD, Lykidis A, Kyrpides NC. 
GenePRIMP: a gene prediction improvement pipe- 
line for prokaryotic genomes. Nat Methods 201 0; 
7:455-457. PubMed 
http://dx.doi.org/10.1038/nmeth.1457 

42. Petersen J. Phylogeny and compatibility: plasmid 
classification in the genomics era. Arch Microbiol 
2011; 193:313-321. PubMed 

43. Petersen J, Brinkmann H, Berger M, BrinkhoffT, 
Pauker O, Pradella S. Origin and evolution of a 
novel DnaA-like plasmid replication type in 
Rhodobacterales. Mol Biol Evol 201 1 ; 28: 1 229- 
1240. PubMed 

http://dx.doi.Org/1 0.1 093/molbev/msq31 0 

44. Petersen J, Brinkmann H, Pradella S. Diversity and 
evolution of repABC type plasmids in 
Rhodobacterales. Environ Microbiol 2009; 
11:2627-2638. PubMed 

http://dx.doi.Org/1 0.1 11 1/j.1 462- 
2920.2009.01987.x 

45. Zielenkiewicz U, Ceglowski P. Mechanisms of 
plasmid stable maintenance with special focus on 
plasmid addiction systems. Acta Biochim Po/2001 ; 
48:1003-1023. PubMed 

46. Petersen J, Frank O, Goker M, Pradella S. 
Extrachromosomal, Extraordinary and Essential - 



The Plasmids of the Roseobacter Clade. Appl 
Microbiol Biotechnol 201 3; 97:2805-281 5. Pub- 
Med http://dx.doi.Org/1 0.1 007/S00253-01 3-4746-8 

47. Cascales E, Christie PJ. The versatile bacterial type 
IV secretion systems. Nat Rev Microbiol 2003; 
1:137-149. PubMed 
http://dx.doi.org/10.1038/nrmicro753 

48. Schwarz S, Hood RD, Mougous JD. What is type 
VI secretion doing in all those bugs? Trends 
Microbiol 201 0; 1 8:531-537. PubMed 
http://dx.doi.Org/1 0.1 01 6/j.tim.201 0.09.001 

49. Chang YJ, Land M, Hauser L, Chertkov O, Larimer 
F, Jeffries CD, Glavina del Rio T, Nolan M, 
Copeland A, Tice H, etal. Non-contiguous fin- 
ished genome sequence and contextual data of the 
filamentous soil bacterium Ktedonobacter 
racemifer type strain (SOSP1-21 7 ). Stand Genomic 
Sci 201 1 ; 5:97-1 1 1 . PubMed 
http://dx.doi.Org/1 0.4056/sigs.21 1 4901 

50. Chan KL, New D, Ghandhi S, Wong F, Lam CMC, 
Wong JTY. Transcript levels of the eukaryotic trans- 
lation initiation factor 5A gene peak at early G1 
phase of the cell cycle in the dinoflagellate 
Crypthecodinium cohnii. Appl Environ Microbiol 
2002; 68:2278-2284. PubMed 
http://dx.doi.Org/1 0.1 1 28/AEM. 68.5.2278- 
2284.2002 

51 . Altschul SF, Gish W, Miller W, Myers EW, Lipman 
DJ. Basic local alignment search tool. / Mol Biol 
1990; 215:403-410. PubMed 

52. Altschul SF, Wootton JC, Gertz EM, Agarwala R, 
Morgulis A, Schaffer AA, Yu YK. Protein database 
searches using compositionally adjusted substitu- 
tion matrices. FEBS J 2005; 272:51 01 -51 09. Pub- 
Med http://dx.doi.Org/1 0.1 1 1 1/j.1 742- 
4658.2005.04945.x 

53. Auch AF, Von Jan M, Klenk HP, Goker M. Digital 
DNA-DNA hybridization for microbial species de- 
lineation by means of genome-to-genome se- 
quence comparison. Stand Genomic Sci 201 0; 
2:117-134. PubMed 
http://dx.doi.org/10.4056/sigs.531120 

54. Auch AF, Klenk HP, Goker M. Standard operating 
procedure for calculating genome-to-genome dis- 
tances based on high-scoring segment pairs. Stand 
Genomic Sci 2010; 2:142-148. PubMed 
http://dx.doi.Org/1 0.4056/sigs.541 628 

55. Meier-Kolthoff JP, Auch AF, Klenk HP, Goker M. 
Genome sequence-based species delimitation with 
confidence intervals and improved distance func- 
tions. BMC Bioinformatics 2013; 14:60. PubMed 
http://dx.doi.Org/1 0.1 1 86/1 471 -21 05-1 4-60 



http://standardsingenomics.org 



141 



